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Abstract. We consider a quantum Hall system with an antidot acting as an energy dependent 
scatterer. In the purely charge case, we find deviations from the Wiedemann-Franz law that take 
place in the nonlinear regime of transport. We also discuss Peltier effects beyond linear response 
and describe both effects using magnetic-field asymmetric transport coefficients. For the spin 
case such as that arising along the helical edge states of a two-dimensional topological insulator, 
we investigate the generation of spin currents as a result of applied voltage and temperature 
differences in samples attached to ferromagnetic leads. We find that in the parallel configuration 
the spin current can be tuned with the leads’ polarization even in the linear regime of transport. 
In contrast, for antiparallel magnetizations the spin currents has a strict nonlinear dependence 
on the applied fields. 


1. Introduction 

Under the application of strong fields, electron transport becomes nonlinear and new effects 
arise: rectification [U 0 El HI 0 El [7], magnetic-field asymmetries [0 0 [TUI El El El El 
El El El El El El El El El E], and the generation of higher harmonics [20 El 127] . 
Quantum coherent conductors are excellent platforms to observe and study nonlinearities since 
small voltage biases cause a sizeable effect over rather short lengths. In the ballistic regime of 
transport, current response is determined by the transmission probability t(E). When t depends 
weakly on the carrier’s energy E, current is always a linear function of V independently of the 
background temperature T. Quite generally, however, transmission across nanostructures shows 
sharp features arising from strongly energy dependent t(E). Furthermore, in the nonlinear 
regime of transport, the transmission becomes a function of the applied voltage V because the 
applied held modifies the potential landscape which in turn affects the scattering properties of the 
sample [28. 29]. The self-consistent procedure must then include electron-electron interactions 
that restore current conservation and gauge invariance beyond linear response. In fact, it is the 
nontrivial dependence of the screening potential on voltage and magnetic fields that explains 
both rectification effects and magnetic-asymmetries, respectively. 

Electron motion can also be induced with the application of external thermal gradients. The 
two main thermoelectric effects are the Seebeck and Peltier effects. The former leads to the 
creation of a voltage drop in response to a temperature difference 0 in the open circuit case. 
The latter is based on the fact that electrons carry energy in addition to charge, and a heat 
current then flows in the presence of an electric current. Both phenomena have been observed 
in mesoscopic systems [m El E], and the agreement with the scattering approach [33] is 


remarkable. In these nanodevices, the thermodynamic efficiency can be tuned with an external 
magnetic field [SUMIETUMUMIIIQ]. 

Beyond linear response, one must also take into account that the transmission is not only 
a function of E and V but also depends on 6: t(E,V,9) [H]. As a consequence, the injected 
charge that builds up in the vicinity of the sample is determined from both particle and entropic 
injectivities m sa ss]. Upper bounds to the performance of heat engines and coolers are thus 
to be carefully evaluated [TF Surprisingly, the thermovoltage can vanish for a finite value 
of 9, as recently demonstrated in experiments with nanowire quantum dots gB]. Furthermore, 
reciprocal relations are shown to break down for sufficiently high heating currents applied to 
multiterminal setups m- Interestingly, departures from the Onsager-Kelvin relations differ 
from auto- or cross-terminal measurements gnssi- 

In this work, we discuss thermoelectric effects beyond linear response for quantum Hall bars 
with an inserted antidot. This system allows for an investigation of the underlying symmetries 
while providing, at the same time, a strong interaction and energy dependent scattering g9ll50] . 
We find large deviations of the Wiedemann-Franz law and interesting nonlinear behavior of the 
Peltier effect. Additionally, we discuss the spin transport in quantum spin Hall systems m 
coupled to ferromagnetic leads. Whereas the normal case already leads to the generation of 
pure spin currents [52], we here find a competition with the polarized currents injected from the 
magnetic contacts. 

2. Theoretical formalism 

When a mesoscopic conductor is coupled to multiple terminals a, (5,, each terminal can be 
characterized by voltage bias eV a = p a — Ep (/ i a is the electrochemical potential and Ep is the 
Fermi energy) and also by the temperature gradient 9 a = T a — T (T a and T are the reservoir 
and the background temperature, respectively). In the presence of these two driving fields, the 
electronic and heat transport is entirely described by the scattering matrix s a p = s a p(E, eU), 
which is a function of the carrier energy E and the potential landscape U inside the conductor. 
The potential U = U(f, {U 7 }, {# 7 }) is, in general, a function of the position r and the set of 
applied voltages {Vy} and thermal gradients {0 7 }. The charge and heat currents at lead a are 
respectively given by mm 
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where A a p = Tr [8 a p — s' a pS a p\ and fp(E) = (1 + exp [{E — up)/ksTp])- 1 is the Fermi distribution 
function in the reservoir fd. In the weakly nonlinear regime of transport, one can expand these 
currents around the equilibrium state, i.e., /j, a = Ep and T a = T, up to second order of V a and 
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The linear conductance coefficients are given by [33] G a p = (2 e 2 /h) f dEA a p(E)(—d E f) ~ 
(2 e 2 /h)A a p(E F ), L a p = (2 e/hT) f dE(E-E F )A a p(E)(-d E f) « (2eir 2 k 2 B T/3h)d E A a p(E)\ E=EF , 
R a /3 = (2 e/h) f dE(E - E F )A a p(E)(-d E f) « (2eir 2 k 2 B T 2 /3h)d E A a p(E)\ E=EF , and K a p = 


(2 /h) f dE ^ E ~ Ep ^ A a p(E)(—dEf ) ~ {2n 2 k^T/3h)A a p(EF), where we have resorted to the Som- 
merfeld expansion for the approximate expressions. We emphasize here that the linear coeffi¬ 
cients G a p, L af 3 , R a p, and K a g are equilibrium quantities and hence are independent of the 
screening potential U. In contrast, the nonlinear coefficients G a p 7 , L af j 7 , M Q/ g 7 , R a p 7 , Kap-yi 
and H a fj-y manifestly depend on U in response to the applied voltage and temperature biases. 
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Leading order interaction effects can be incorporated into the nonequilibrium potential U via 

characteristic potentials (CPs) u a = ( dU/dV a ) eq and z a = (dU/dd a ) eq : 


U — U e q + ^ ( U a V a -p ^ [ Z a 9 c 


(5) 


Note that CPs u a and z a relate the variation of U to voltage and temperature shifts, 
respectively. We determine U self-consistently by considering the net charge of the system 
q = Qbaie + Qsci- The bare charge gWe = e Ea(^« e hc* + D e a 6 a ) has two contributions due 
to voltage bias and temperature shift, described by the particle injectivity [28] [29j v^{E) = 

(2vr i)- 1 E/3 Tr [4 a^mr] and the entropic injectivity [Tlj K( E ) = ( 2?r *) _1 E /3 Tr 
with Da ,e = — f dEva e (E)dEf ■ The screening charge q SCT building up inside the sample can be 
obtained from the response of the screening potential, A U = U — U eq , away from the equilibrium 
state U eq . One can write q SCT = e 2 IIAf7 where II is the Lindhard function which becomes 
II = f dED(E)dEf in the long wavelength approximation, with D{E) being the density of states. 
Then, the set of equations for the CPs is closed via the Poisson’s equation V 2 A U = —47rg. 

Finally, for a practical calculation, one can resort to the WKB approximation valid in the long 
wavelength limit and make the replacement S/SU -A —ed/dE, i.e., 9 gy^ = ~ —eu 7 d Q^ 

and 9 qq P = ~ ~ ez j~§^- Tor a two-terminal setup which we consider below, the matrix 

elements of A a p are given by An = A 22 = —A 12 = — A 21 = t(E), where t(E) is the transmission 
probability through the system. 


2.1. Magnetic-field asymmetry 

At equilibrium, the screening potential U = U eq is symmetric with respect to the reversal of an 
applied magnetic field, i.e., U eq (B) = U eq (—B), due to the fundamental microscopic reversibility 
principle. In Ref. [33j, magnetic-field symmetry of linear thermoelectric and heat transport 
has been shown in the framework of scattering theory. When the system is driven into the 
nonequilibrium regime, however, this magnetic-field symmetry can be broken since the CPs in 
Eq. dSJ) are in general magnetic-field asymmetric, i.e., u a (B) 7 ^ u a (—B) and z a (B) z a {—B). In 
Ref. [S], the nonlinear conductance G a g 7 in the isothermal case has been shown to be magnetic- 
field asymmetric since u a , the voltage response of U, is not an even function of B. Here, we 
generalize this result to all thermoelectric and heat transport coefficients. 

We quantify the magnetic-field asymmetry in the nonlinear transport regime by defining the 
symmetry (E) and the asymmetry (A) parameters [38] for G , L, R , and K coefficients in Eqs. ([2]) 
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where X a p indicates the differential transport coefficients Gap, C a p, 7 Z a p, and K a p defined by 
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and ^“ ear refers to the linear terms G a p, L a p, R a p, and R a p- We here consider either an 
isothermal, i.e., { 0 } = 0 , or an isoelectric case, i.e., {Id} = 0 , hence the terms XI a pp and H a p 1 in 
Eqs. m and (HD are not considered. Notice that X a g contains both linear and nonlinear terms 
and in the linear response regime it satisfies ^^p pa = Xyp p a = 1 > due to the microscopic 
reversibility X^ ear (B) = A’g 1 ^ ear (— B). Hence, a deviation from 1 indicates the magnetic- 
field asymmetry in the nonlinear regime. Moreover, these definitions have direct relevance to 
experiments 03 and also are related to the efficiency of the thermoelectric power generation or 

the refrigeration [331132 Eg]. 


3. Quantum Hall bar wih an antidot 

As a model system for the general formalism described above, we consider a conductor in the 
quantum Hall regime in a two-terminal setup, see Fig. [I] We suppose that B is strong enough so 
that only the lowest Landau level is occupied and consider the reversal of its direction B —>• —B. 
A gate-controllable antidot [19, hO] can connect two counter-propagating edge states, which we 
regard as a quantum impurity with a Breit-Wigner resonance at Ed + eUd(B ) where Up(B ) is 
the potential shift at the antidot in the presence of magnetic field B. The upper and lower edge 
states are tunnel-coupled to the antidot level via coupling strengths Ti and 17 , respectively. 
If the direction of the magnetic field is reversed, the resonant level at the antidot is located 
at Ed + eUd{—B). It should be noted that in general Ud(B ) pP Ud(—B) beyond the linear 
response regime [8]. The H-asymmetry can appear either via scattering asymmetry, Ti ^ T 2 , or 
electrical asymmetry if the Coulomb interaction between the upper edge and the antidot charges 
is stronger (or weaker) than that between the lower edge and antidot. 

In Fig. [33 we divide the conductor potential into five regions f \ with i = 1,..., 5, with H 5 = 
indicating the antidot region. We assume that the potential Ui in each region is constant and the 
Coulomb interactions between different regions are described by a capacitance matrix Cij [28]. 
capturing the essential physics of our interest [§U2B1129j. For definiteness, we also assume the 
equal density of states for all regions, i.e., Di = D and the symmetric injectivities between two 
terminals, i.e., H?’ 6 = D p,e and n* = n. Finally, we self-consistently determine the internal 
potential and therefore the CPs by solving the Poisson’s equation explained previously. 

We here consider the two cases: (i) the scattering asymmetry, i.e., Ci = C,T 1 = (1 + ??)r/2, 
r 2 = (1 — r/)T/2, (ii) the electrical asymmetry, i.e., Ti = p 2 , C\ = C 2 = (1 + £)(7, 
C 3 = C 4 = (1 — £)C. In the above two cases, the inherent asymmetry is described with a 












Figure 1. A two-terminal quantum Hall bar-antidot system, where voltages Vi, V 2 and 
temperature gradients 6\, do are applied. An antidot (H 5 ) is coupled to the edge states with the 
hybridization widths Pi and T 2 and the Coulomb interactions are described by capacitances C \, 

d, c 3 , c 4 . 


parameter 77 or £, respectively. We evaluate the CPs as 
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where r = 1 — t = r iP 2 /|A| 2 is the Breit-Wigner reflection ( t : transmission) through the antidot 
at equilibrium, and A = Ep — Ed — et/^ q + *T/2. 

In Eq. (HD, one can easily prove the fundamental sum rule u\{B) + u 2 (B) = u\(—B) + 
u 2 {—B) = 1 due to the gauge invariance. Importantly, the CPs are generally magnetic-field 
asymmetric, i.e., u a (B) / u a (—B) and z a (B ) / z a (—B). We also point out the properties 
ui(±B) = u 2 (^fB) and z\ (EB) = z 2 {^B) in Eqs. (fill and (HD, due to the chiral nature of the 
quantum Hall system. 

The symmetry (E) and the asymmetry(A) parameters defined in Eq. Q can be evaluated 
straightforwardly with the CPs in Eqs. (HD and ( 1121 ) . The general expressions of these 
parameters for a two-terminal quantum conductor can be found in Ref. [15], to which we also 
refer the readers for the detailed analysis of our quantum Hall system. We here briefly summarize 
some of the interesting results. 

Firstly, we find all the off-diagonal asymmetry parameters as well as the symmetry 

parameter n for electric conductance always exhibit H-symmetry even in the nonlinear 
regime: 
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This originates from the chiral property u\{EB) = u 2 (^B) and Z\{EB) = z 2 (^B). In addition 
to this, the gauge invariance condition u a (±B) = 1 gives a contribution for the derivation 
of S| j , j = 1 because u\(B) + u\{—B) = u\(B) + u 2 (B ) = 1 where the latter equality comes 
from the gauge invariance (see Ref. [48] for the details). 























Secondly, we find that the symmetry parameters for £ (thermoelectric) and /C (thermal) 
coefficients depend on the lead indices: 


Sfi,u = 1 + cf (20/T), Sf 2>21 = 1 + cf (2 9/T), 

Sfi,u = 1 + cf(20/T), Sf 2 , 21 = 1 + c%(29/T). 


(15a) 

(15b) 


This different tendency between diagonal and off-diagonal parameters have a relevance to recent 
experiments 023- Intriguingly, we note that this is a high-temperature effect since in the limit 
ksT —> 0 we find Sfj n = Ef 2 2 \ = ^n,n = E)f 22 1 = 1 + 2 6/T, irrespective of the system 
parameters (see Ref. [32] for the details). 

In this quantum Hall system, eight parameters n and Ay 2 21 for all X = Q, £, 1Z, /C, 
are independent of the scattering asymmetry^) and the electrical asymmetry(£) factors due to 
the chiral nature, including five R-symmetric parameters in Eq. (HD- However, the off-diagonal 
elements Ef 2 21 an d Ef 2 21 are dependent on r/ or £ Consequently, the distinction between the 
diagonal and the off-diagonal elements disappears for ?y = £ = 0. Thus, an asymmetry inherent 
in the system is crucial to observe this distinction, which is consistent with an asymmetric 
scattering in a recent experiment 03- 

Even with a nonzero rj or £, our results show that we can tune the antidot level to make 
cf = c 2 = 1 in Eq. (1 15 all and recover the universality of the thermoelectric coefficients, i.e., 
Sfj j | = Ef 2 21 = 1 + 2 6/T. But, this is not the case for the heat current counterparts Ef) n 
and Ef 2 , 2 i- 

In addition to the above parameters, one can also write E 22 2 i = 1 — (P{eV/T), A 21n = 
1 + c g (eV/r), Af l u = l + c%(20/T), Af x u = 1 + c%{eV/T), and Hf in = l + c$(2 9/T), where 
(P = cf = = Pf = 0 if ?? = £ = 0. Hence R-asymmetry of these parameters is due only to 

the underlying asymmetry in the system. We note here that Ey 2 21 and Afj n are described by 
a single constant P. 

Finally, we write the electrothermal symmetry parameters E(f n = 1 + c^feV/T) and 
Ef 2 21 = 1 — cp(eV/T) in heat current measurements. Here, E^f 21 is in general a function 
of 77 or £, whereas E(f n is not. 

3.1. Nonlinear Peltier effect 

When an electrical current (I) flows through a conductor in a isothermal configuration a heat 
current (J) is generated. The linear Peltier coefficient is defined as the ratio between the heat 
and electrical currents, 

n 0 = , (16) 

In Ref. [13] the n coefficient was generalized to the nonlinear regime. By keeping only the 
leading order corrections in powers of the electrical current one has 
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The conversion factor for electric currents into heat flow in the nonlinear regime is given by 
the relative strength of the nonlinear conductances, Rm, and Gm to the linear ones (R 11 , and 
Gu). In terms of the symmetries of the nonlinear coefficients we can rewrite Eq. (1171) as 
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In Fig. [2] we plot dn normalized to a given value of the electric current /. We observe large 











Figure 2. Plots of <511 versus antidot level £q/T for several background temperatures kpT. We 
used kpO/T = eV/T = 0.1 and rj = 0.1. 


deviations of the Peltier effect arising around resonance (i.e., when eo ~ Ep). Otherwise, the 
nonlinear Peltier coefficient attains its linear value. Besides, we find that deviations occurring 
for eo ~ Ep are almost independent of the background temperature. Only when eo lies off 
resonance, the largest deviations take place for the highest temperature. 


3.2. Wiedemann-Franz law 

We now explore the degree of fulfilment of the Wiedemann-Franz law beyond linear response. 
For a metallic system and at very low temperatures the ratio between the thermal and electric 
conductances normalized to the background temperatue T equals a universal value, 


A 0 = 


A'n 
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where Ao = yy(^r) 2 is the Lorentz number. In the nonlinear regime, Ref. [53] suggests a 
similar quantity as the ratio between the heat flow (normalized to the temperature shift) in the 
isoelectrical configuration and the electrical current in the isothermal case (normalized to the 
applied voltage), 
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The normalized deviation of the Wiedemann-Franz law from the Lorentz number can be 
expressed in terms of the nonlinear transport coefficients, 
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which, in turn, can be recast using the symmetry and asymmetry coefficients defined above, 
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Figure 3. Plots of (A — Ao)/Ao versus antidot level Eq/T for several temperatures ksO. We use 
ksT/T = 0.1, eV /T = 0.05 and g = 0.5. 


We find that departures become stronger for level resonances around the Fermi energy within 
a energy scale given by T. This is so because the T determines the energy variation of the 
transmission function. In addition, the deviations increase for higher 6 since the heat transport 
is more linear with increasing 6. 

4 . Thermoelectric transport in quantum spin Hall insulators 

4-1. Spin-generalized screening potential 

In the previous section, a quantum Hall antidot in a two-terminal setup has been investigated, 
where we determine the CPs by solving the Poisson’s equation. Here, we consider a simple 
extension to the quantum spin Hall (QSH) system, see Fig.[H In order to deal with the spintronic 
case with the same methodology, we generalize the scattering formalism in a spin-dependent 
manner. The potential U a = U(f, {Ey}, {# 7 }, a) is now a function of the spin index a =t,i 
as well. This cr-dependence is crucial in our QSH system due to the underlying helicity, i.e., 
the spin-momentum correlation. Indeed, the property U ^ ^ Ui through the antidot filter is the 
operational principle for the spin-polarized currents even with normal metallic contacts. 

We generalize Eqs. © and (O into the spin-resolved form |52| 

= \ dEA°p(E,eU)fp(E), 

P J 

X = T E / dE ^ E - ^)A a a p(E,eU)fg(E), 

(3 J 

for which we divide 2 A a p in Eqs. (JT]) and (J2J) into A^ = A a g(U^) and A^g = A a g(Ui), 
explicitly incorporating the spin-dependent screening effect. In a two-terminal setup, we have 
Ah = A %2 = —Af 2 = —A -21 = t a (E) where t a (E) is the spin-resolved transmission probability. 
We disregard the spin-flip scattering, hence the current conservation condition is satisfied for 
each spin component as Yla^a = 0 and + d^Va) = 0- One can define the direction 
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Figure 4. Schematics of quantum spin Hall bar antidot setup. Spin-generalized interactions 
are described by spin-dependent capacitance Ci St da, with i = 1,2 being the edge label, s = ± 
the helicity, d standing for the dot, and a =t,4- the electronic spin. Spin-dependent couplings 
between the helical edges and the single level antidot are denoted with Tj s . 


of spin-resolved currents I a = /f = —If and = —J% — I a (Vi ~ V 2 ), and hence the 

spin-polarized currents as well: 


h = i\-h 

Js = J\ — J\ 

Total charge and heat fluxes are given by I C = I^ + 1^ and J c = respectively. 

The screening potential U = Yh a U a in Eq. d5]) can now be generalized as 
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where u aa = (dU a / dV a ) e(i and z aa = (dU a /d6 a ) eq are spin-dependent CPs which relate the 
variation of spin-resolved potential U a to voltage and temperature shifts at each terminal. 

As described earlier, the self-consistent determination of U can be accomplished via Poisson’s 
equation V 2 AC = —47rg, with Af7 = U — f7 eq = A U a and 


q = ^2q a = eY^ \D p a (cr)eV Q + D e a (a)0 a 

cr a,a 


+ e 2 ^n CT A U a . 

<J 


(28) 


In our model, cr-dependences of Da' e (cr) and appear for unequal spin populations arising, 
e.g., from ferromagnetic leads. In order to understand this, note that the first two terms in a 
square bracket in Eq. (1281) are the contributions from the lead injection of charges by means of 
voltage and thermal driving, respectively. Only when the spin population is unequal (p / 0, 
p: polarization) via using ferromagnetic contacts, one expects Da' e ( t) 7 ^ -Da’ e ( 4 -) and n>|- + n ; , 
where the latter comes directly from D + 7 ^ D±. For normal metallic contacts, the only term 
in Eq. (1281) which can give rise to a spin imbalance inside the system is the screening potential 
A U a . 

The solution procedure is quite analogous to the quantum Hall case (see Ref. [52] for details). 
It should be noted that for QSH case the edge-antidot couplings Ti s and T 2 S in general depend on 
the helicity s = ± corresponding to spin channels t(+) and j,(—), when coupled to spin-polarized 
ferromagnetic contacts. This is also the case for the reflection and transmission probabilities 
since r a = 1 — t a = r i s r 2s /|As | 2 , where A s = E F — E d + iT s /2 with T s = Ti s + r 2s . But, 
for normal contacts with p = 0, there is no spin imbalance at the edge states, which leads to 




















= ^ 4 . antidot scattering. In this case, the linear conductances are spin-independent and 
the spin-polarization appears only in the nonlinear transport regime. In contrast, for general 
case with p / 0, unequal spin density leads to t^(p) ^ t\.{p) giving rise to spin-polarized electric 
and heat currents already in the linear regime (see Ref. [52J for details). However, even for p / 0 
case, if the two contacts are magnetized in an antiparallel configuration, one finds tf(p) = t±(p) 
and hence the linear spin-polarization vanishes as in normal contacts, due to the helical nature 
of the QSH system. 

The effective Poisson’s equation reads 

Qis = e ’y^ J (Djs.a e V°‘ + + e 2 Hi s AUis = Ci St da{AUj s — A Uda)- (29) 

Ot <J 

Note here that the charge with spin a =t(i) at the antidot is supplied from the edge with helicity 
s = +(—) since in our model we neglect the spin-flip scattering, by which one can maximize spin- 
polarization effects. The density of states for all regions are given by Di S = D s = (1 + sp)D/ 2, 
and the symmetric injectivities from the two terminals give D?’ e Q = Dg’ e = (1 + sp)D p,e / 2 and 
n is = n s = (1 + sp)U/2. 

4-2. Normal contacts 

Here, we focus on the result for normal metallic contacts (p = 0) since this is the most attractive 
case of an all-electrical setup, in which our main physics of interest resides. For p = 0, it 
should be mentioned that the spin-channel density of states does not depend on the helicity, i.e., 
D + = D _ = D/2 with D = D + + D _. 

We consider the scattering asymmetric case, i.e., equal interaction strength Ci St da = C ) 5 = 
C s = C/2 with C = C + + C- , but asymmetric hybridizations = (l + r/)T/4, T 2 S = (1 —r/)T/4 
with T = T + + T_ (r s = + T 2 S = r/2). The coupling asymmetry is thus quantified with 

a nonzero p = (Ti — T 2 )/T where T* = ^ S P; S . In an experiment, this can be the general case 
where the antidot is closer to one of the edge states. Also, this asymmetry can be generated by 
tuning the width and the height of the tunnel barriers between the antidot level and the edge 
channels. One can find A Uda = ui cA'i + U 2 a ^2 + zio-#i + " 2 <t# 2 5 where the corresponding CPs 
are given by 
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Here, c sc = [2 — 2C/e 2 n] _1 = C^/2C with 1 /C ]x = 1/C + 1/e 2 D. Importantly, we can expect 
electronic transport to be spin polarized for asymmetric couplings due to the spin-dependent 
screening. Intriguingly, our result is a purely interaction-driven effect and this effect disappears 
in the noninteracting limit C —> oo. 

Up to the second order expansion of V\ = V and 9\ = 6 (V 2 = O 2 = 0), we have the current 
expressions 
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where t = t(Ep), f = dEt(E)\E=E F , and t" = d B t(E)\E=E F - These are central to our analytic 
results. The spin-polarized electronic and heat currents can be generated if r/c sc 0. Note 
also that spin-polarized currents appear in the nonlinear regime only. Isothermal (8 = 0) and 
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Figure 5. (a) I4h versus k B 8/Y and (b) 0 a d versus eV/T, at Ed = 0.25T for several k B T. Used 
parameters are r) = c sc = 0.5 and E B = 0. 


isoelectric (V = 0) cases are treated in Ref. [52]. Here, we focus on the pure spin currents, i.e., 
I c = 0, J c = 0, by means of thermoelectric Seebeck and Peltier effects. 

In open-circuit conditions, a thermovoltage Vth can be generated in response to a temperature 
bias 8 , making I c = 0. Figure [5J a) shows the numerically evaluated set of biases {#,U} which 
satisfies I c (\th,6) = 0, where the slope can be identified as the Seebeck coefficient. Note that 
the thermovoltage acquires a nonlinear component with increasing 8 [31, 36]. One finds the pure 
spin current expression as a function of thermal gradient only 
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up to leading order in 8. Figure [6][a) shows the numerical plots of pure I s versus 8 beyond 
the quadratic regime (a comparison with the analytical results can be found in Ref. [52] 1. We 
observe that the amplitude of I s shows a nonmonotonic behavior with T, providing another way 
to maximize the effect. 

In our setup, also the pure spin heat flows can be generated using electrical means only. 
By adiabatically isolating the sample, thermal bias 8^ is generated in response to the applied 
voltage V, making J C {V,8 a d) = 0, see Fig. 0(b). We find the pure spin heat current expression 
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up to leading order in V. Corresponding plots are shown in Fig. 0(b) as a function of V. 
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Figure 6. (a) Pure I s versus kpO/Y and (b) pure J s versus eV/T, at Ea/T = 0.2 for several 
kpT. Here, r] = c sc = 0.5 and Ep = 0. 


4-3. Polarized contacts 

For polarized contacts with p / 0, an unequal spin injection from the reservoirs with a majority 
spin density (l+p)D/2 and a minority one (l—p)D/2 would generate spin-polarized currents. In 
this sense, the results are to some extent anticipated and hence less interesting than those of the 
unpolarized case. For an equal and antiparallel magnetization, however, it is interesting to note 
that linear conductances are spin-independent giving rise to vanishing linear spin-polarization. 
This will be explained below. 

Although the QSH bar itself preserves the time-reversal invariance, we can also investigate the 
reciprocity relation by considering ferromagnetic contacts with a polarization p, where the role 
of a magnetic-field reversal (B —> —B) can be achieved by simultaneous reversal of polarization 
p —p and spin index a —> a (a = 1 \l) [53j. As we will explain below, in any case the 
Onsager-Casimir relation is always satisfied for linear coefficients, i.e., X^(p) = X^ Q (—p) for 
X = G,L , R, K , which holds even for p = 0. Beyond the above argument, we here discuss 
about the breakdown of this fundamental symmetry relation in nonlinear regime of transport for 
parallel and antiparallel magnetization respectively, and view the normal contacts as a limiting 
case of these two. 


4-3.1. Parallel magnetization For an equal and parallel magnetization between two reservoirs 
with a polarization p, the density of majority spin-component (say, j') is given by = (1 + 
p)Di /2 and the minority one by = (l—p)D^/2 where l = L,R denotes the left or right leads. 
Then, including the helicity s-dependence, we have F^ s = (l+sp)Tn/2 where i = 1, 2 corresponds 
to the upper or lower counterpart for edge-antidot coupling [see Fig. 0 ], This gives polarization- 
and spin-dependent transmission as t a (p,E) = 16 (E — E^) 2 /[16(E — Ej) 2 + (1 + sp) 2 T 2 ] with 
T = Ti + T 2 . Hence, as expected, spin-polarized currents appear already in the linear regime 
since t^(p) — t±{jp) ^ 0. Moreover, one can notice here that t^(p) = t^{—p) because the relation 
sp = s(—p) is always satisfied in the expression of transmission, with s = +(—) corresponding 
to cr =t (|). Thus, we have G*g(p) = Gj a (-p), L^(p) = L^ a (-p), R^p) = and 

(p) = Kp Q (—p) for a,j3 = 1,2 (we refer the reader to Ref. [52] for explicit expressions 












of all coefficients). This is a manifestation of the fundamental symmetry relations near 
equilibrium [33). 

We calculate the screening potential A U a (p) = u f (p)V\ + u 2 (p)V\ + zf {p)0\ + z 2 (p)$ 2 , where 
the corresponding CPs in the presence of (i) scattering and (ii) electrical asymmetries are given 
by 
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Here, t CT (p) = 1 — r a (p) = !6(Ep — Ed) 2 /[16(Ep — Ed) 2 + {1 + sp) 2 T 2 } and s = ± refers to a =t, I- 
We note that ct c (p) = —ct c (—p) and c^ip) = — c ti(~ an d hence u{(p) = u 2 (—p). 
Furthermore, one can show that u[ =l 2 {jp) — u\ =l2 (—p) oc ?/,£ and u^ =12 (p) ~ u\ =l2 (p) = 
±rj[ct c (p) - cic(p)] or ±f [cjj(p) — c^i(p)]. Importantly, the former characterizes the breakdown of 
Onsager-Casimir symmetry in nonlinear regime [481 15]. e.g., G iu(p)-G\u(~P) oc V,€, w bile the 
latter does the nonlinear contribution of the spin-polarized currents, e.g., Gj n (p) — G\ { j (p) = 
g(r],^;p), with g(r],£ m ,p) being a function of 77 or £ and p. For p = 0, these two relations merge 

f I I 

into a single one u\ — uj = 2r/c sc , 2£c e i, leading to G{ 11 — G\ lx oc 77 , £. Previously, this has 
been shown to generate spin-polarized currents in the nonlinear transport regime with normal 
metallic contacts. 
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Figure 7. Plots of I s /I c versus (a) eV/T at kpT/T = 0.01 and (b) kpO/T at kpT/Y = 0.5, for 
several (parallel) polarization p, with 77 = c sc = 0.5, Ed/T = 0.25, and Ep = 0. With p ^ 0, 
spin-polarization occurs even at zero biases, i.e., V = 0 or 9 = 0. 


















Figure 8 . Plots of J S IJ C versus (a) eV /F at kpT/T = 0.01 and (b) kpO/T at kpT/T = 0.05, 
for several (parallel) polarization p, with 77 = c sc = 0.5, E^/T = 0.3, and Ep = 0. With p / 0, 
spin-polarization occurs even at zero biases, i.e., V = 0 or 9 = 0. 


4-3.2. Antiparallel magnetization For an equal but antiparallel magnetization between two 
reservoirs, however, the upper edge is occupied by right-moving and left-moving ^-electrons 
both with an equal portion (1 +p)/2 of the total density of states in each lead, because in 
this configuration the majority spin component is up (down) for the left (right) reservoir. In 
contrast, the lower helical edge is filled only with a minority portion (1 —p)/2 for both left- 
moving [+(t)] and right-moving [—(4-)] channels. The net effect of this antiparallel configuration 
gives spin-independent coupling and transmission, i.e., Ti s = (1 + p)T\/2, T 2 S = (1 — p)T 2 / 2 , 
and t a {p , E) = [16(E — Ed) 2 +p 2 T 2 ]/[l6(E — Ed) 2 + T 2 }. We then have t^(p) = t±(p), hence linear 
conductance coefficients make no contribution to spin-polarized currents, i.e., G^(p) = (p ), 

L ip(p) = L t/ 3 (P)> R !p( p) = ^(P)> and K lp(p) = K tp(p)- This vanishing linear order 
contribution even with a nonzero p but with a specific antiparallel configuration is due to 
the helical nature of the QSH system. Nevertheless, the symmetry relation still holds, i.e., 

Mp) = h(-p)< llence G lp(p) = L Up) = L l !„(-?). Kfl(p) = and 

Kl„(p) = K^(-p) 

We analogously evaluate A U a (p) = u" (p) Vj + u 2 (p)V‘i + zf (jp)0\ + z%(p)02, where 
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Here, t(p) = 1 — r(p) = [16 (Ep — Ed) 2 + p 2 r 2 ]/[16(Ep — Ed) 2 + T 2 ] and s = ± corresponds to 
cr =t,|- 

















It is easy to show that ct c (p,p) = —ci c (r?,p), c^(£,p) = -c^fop), ct c (p = 0 ,p) = ct c (p = 
0 ,-p), 4(£ = 0 ,p) = = 0,-p), < c (r/ = 0,p) oc sp, c&(£,p) oc sp, c° c {p,p = 0) oc rjs, and 

u°ip) =u%[p). Thus, we have, e.g., u\{p)-u\(-p) = cL(p,p)~cl c (p, -p) = c£ c (r],p) + ct c (p, -p) 
as well as u[(p) - u\(p) = ct c (p,p) - cj c (p,p) = 2ct c (rj,p). 

The relation such as u\(p) — u\(—p) = ct c (p,p) + ct c (p,— p) indicates the Onsager-Casimir 
symmetry breaking in nonlinear regime, which vanishes for p = 0 but survives for p = 0 , 
p 7 ^ 0 with the latter corresponding to the normal metal leads. The nonlinear spin-polarization 
term, e.g., u\(p) — u\{p) = 2ct c {p,p) remains finite provided that either of p or p is nonzero. 
For p = 0, as in the parallel configuration, these two relations merge into a single one, e.g., 
u[ — u\ = 2pc SCl 2^c e \. Only when p ,£ = p = 0, all spin-polarized currents vanish, which again 
explains the spin-filter effect in an unpolarized case [Eqs. (1321) and (133]) ] . as a limiting case p — > 0 
of the antiparallel configuration. 

From the symmetry arguments in this section, one can notice that the Onsager-Casimir 
symmetry breakdown in nonlinear regime for p —> 0 with a nonzero asymmetry factor p or £ 
suggests the underlying principle for the observed spin-polarization for QSH antidot system. 

5. Conclusions 

We have analysed the nonlinear transport of quantum Hall setups subjected to the voltage and 
temperature biases. Our theory is based on the scattering transport formalism that incorporates 
electron-electron interaction within a mean-field description. The potential landscape of the 
nanostructure depends on the injected charges due to both, the voltage and temperature 
gradients. Using this theory we find the symmetry relations of the nonlinear transport 
coefficients for a quantum Hall bar with an inserted antidot. We find large deviations for the 
Wiedemann-Franz law and Peltier coefficient from their linear values. We also have reported 
results for the spin Hall bar including the antidot system in a spin-dependent capacitative 
model. We have found pure spin current generation for which the spin energy and spin 
current values depend on the magnitude of the lead polarization. Our results are relevant 
in view of recent works that emphasize nonlinear properties of thermopower in nanostructures 
[54l 15511561157115811591 [60l I 6 T 11621 1631 [Ml US]- In general, the nonlinear regime of transport has 
not been explored to a large extent. Therefore, we expect the coming years to be full of exciting 
discoveries and developments in this held. 
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